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Abstract 

In a one-pion plus one-kaon exchange model, we calculate the kaon production 
cross sections in nucleon-nucleon, nucleon-delta and delta-delta interactions. 
We find that this model describes reasonably well the experimental data on 
kaon production in the proton-proton interaction. Near the kaon production 
threshold, the cross section obtained from this model is smaller than that 
from the linear parameterization of Randrup and Ko. For kaon production 
cross sections from the nucleon-delta and delta-delta interactions, the cross 
sections are singular in free space, so we calculate them in a nuclear medium 
by including the (complex) pion self-energy. The results are compared with 
the scaling ansatz of Randrup and Ko. The theoretical cross sections are then 
used in a transport model to study kaon production from Au+Au collisions 
at 1 GeV/nucleon. 
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I. INTRODUCTION 



Kaon production in heavy-ion collisions has been extensively studied both experimentally 
and theoretically |j^ 15,?, 17-21 1 at various incident energies. Its interest stems from the 
fact that because of strangeness conservation, kaons will not be absorbed in nuclear medium 
once they are produced in the collisions. They are thus expected to carry useful information 
about the early stage of heavy-ion collisions. Indeed, at incident energies below the kaon 
production threshold in the nucleon-nucleon interaction in free space (~ 1.58 GeV), the kaon 
yield has been shown in transport models to be sensitive to both the nuclear equation of 
state at high densities ||7|, p!0| , [Tll , [T5H l7|l and the kaon scalar potential in nuclear medium [[13 



The latter is related to the partial restoration of chiral symmetry |]22| , p3| . From comparing 
theoretical results with experimental data measured in Au+Au collisions at 1 GeV/nucleon, 
it has been concluded that the attractive kaon scalar potential in dense nuclear matter is 
appreciable and the nuclear equation of state at high densities is soft. 

For heavy ion collisions at incident energies of around 1 GeV/nucleon, the colliding 
system consists mainly of nucleons, deltas and pions. Kaons can be produced from baryon- 
baryon (nucleon-nucleon, nucleon-delta, and delta-delta) and pion-baryon (pion-nucleon and 
pion-delta) interactions. Kaon production cross sections in these elementary processes are 
thus needed in transport models in order to evaluate the kaon yield in heavy-ion collisions. 
The most straightforward way is to use experimentally measured cross sections. There exist 



some experimental data for kaon production cross sections from the pion-proton and 
proton-proton interactions. For the proton-proton interaction, there are no data near the 
the kaon production threshold of about 1.58 GeV, which are the energies at which kaons are 
mostly produced in heavy ion collisions at around 1 GeV/nucleon. The extrapolation from 
the available experimental data at high energies to near the threshold thus depends on the 
parameterization used for the cross sections. There are two popular parameterizations of 
the experimental data; the linear parameterization of Randrup and Ko |^ and the quartic 
parameterization of Zwermann B. Although both fit the experimental data at the lowest 
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available incident energy, they differ considerably near the kaon production threshold. 

Furthermore, the kaon production cross sections in nucleon-delta and delta-delta in- 
teractions are not available experimentally. In the work of Randrup and Ko [^], these 
cross sections are related to that from the nucleon-nucleon interaction mainly through argu- 
ments based on isospin symmetry. This leads to a simple scaling ansatz: <Jnn^nyk{^/s) = 
(4/3)(jArA^Aryii-(-\/s) = '2'0'aa->nyk{^/s)- No detailed theoretical study has been carried out 
to examine its accuracy. All existing calculations based on transport models have shown 
that deltas play a significant role in subthreshold kaon production due to their larger masses 
than the nucleon mass [0,|l^,|ll|,|l^ . In these calculations, the scaling ansatz of Ref. 



have usually been used. At higher incident energies in AGS experiments, deltas have also 
been found to be important for kaon production. In microscopic models such as the ARC 
19| and the RQMD |2^, the kaon production cross sections from nucleon-delta and delta- 



delta interactions are assumed, however, to be the same as that from the nucleon-nucleon 
interaction. 

The above discussions thus point to the need for a theoretical model in evaluating the 
elementary kaon production cross sections in hadron-hadron interactions. In Ref. PSI, a 



resonance model has been used to study kaon production from the pion-nucleon interaction 
and has been shown to reproduce rather well the available experimental data. The resonance 
model has also been extended to kaon production from the pion-delta interaction. For kaon 
production in the nucleon-nucleon interaction, a one-pion exchange model has been used 
in Ref. p6[ by including also medium modifications of the pion propagation in the delta- 



hole model. The contribution of one-kaon exchange to kaon production in proton-proton 



interaction has been studied in Ref. ^7j, and has been found to be as important as the 
one-pion exchange contribution if off-shell effects are included. No detailed calculation has 
been carried out for kaon production from the nucleon-delta and delta-delta interactions. 

In this paper, we shall calculate kaon production cross sections in nucleon-nucleon, 
nucleon-delta and delta-delta interactions based on a one-pion plus one-kaon exchange 
model. In Section 2, we present the formalism for these cross sections. The results are 



presented in Sections 3 and compared with available experimental data and existing param- 
eterizations. In particular, we discuss kaon production cross sections from the nucleon-delta 
and delta-delta interactions and compare them with the scaling ansatz of Ref. 0. In section 
4, we use the theoretical elementary kaon production cross sections in a transport model 
to study kaon production from Au+Au collisions at 1 GeV/nucleon. Finally, a summary is 
given in Section 5. 



II. KAON PRODUCTION IN BARYON-BARYON INTERACTIONS: 

FORMALISM 

In this section, we present the formalism for kaon production cross sections in nucleon- 
nucleon, nucleon-delta and delta-delta interactions. Only final states without a delta will be 
considered as the threshold energies for final states with a delta are much higher than the 
average energy of two interacting baryons in heavy ion collisions at 1 GeV/nucleon. 



A. NN NYK 

The relevant pion and kaon exchange diagrams are shown in Fig. 1, where Y , denoted by 
thick solid lines, represents either a A or a S hyperon. Let us first consider the case in which 
the kaon is produced in association with a A hyperon. To fix the relative sign between the 
pion and kaon exchange amplitudes requires a well defined interacting Lagrangian. Since 
we shall treat the tcNYK and KNNK interactions phenomenologically using their cross 
sections determined from the resonance model |2^, their relative sign can not be fixed. 



Instead choosing the sign to maximize the cross sections as in Ref. [|2^, we shall simply 



neglect the interference between the two amplitudes and fit the experimental cross sections 
by adjusting other parameters in the model. The isospin-averaged cross section is then given 

by mm 



(v^) 
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CTKN^KNiwK)- (1) 



In the above, masses of the nucleon, pion, kaon, and A hyperon are denoted by m^v, /Btt, mx-, 
and mA, respectively. f-nNN and Jkn^ are the pseudovector ttNN and KNA couphng 
constants. The four-momenta of two initial nucleons are denoted by pi and p2, respectively. 
Their total energy in the center-of-mass system is denoted by ^/s, while the magnitude of 
the three-momentum of each nucleon in the center-of-mass frame is p. denotes the total 
energy of the pion-nucleon system in the center-of-mass frame, with {wT^)min — f^K + fn\ 
and {w.,^)max — \fs — mN- Similarly, wr is the total energy of the kaon-nucleon system, with 
{wK)min = t^K + '<^N and {wxjmax = y/s — m\. The magnitude of the three-momentum of 
the exchanged pion, A;^, is related to w^^ by 



and a similar expression for kx, with and in Eq. (2) replaced by uik and wk, 
respectively. 

In Eq. (1), the four-momentum of the exchanged pion is denoted by ^tt- In the center- 
of-mass frame, we have 



where p' is the magnitude of the three-momentum of the final nucleon and is related to ^/s 
and u^TT by 




TT 




{ql)+ = 2m% - 2{m% + pY^\m%+py/^ - 2pp', 



(3) 




(4) 



Similarly, for (q'^)± we have 



= m% + ml- 2{m% + p'f'^ml + p'^'^ + 2pp\ 
{ql)+ ^m% + ml- 2{m% + p^/^ml + p'^/^ - 2pp', 



(5) 
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with p' calculated from Eq. (4) by replacing mjv with mA- 

The pion propagator D^{ql) and the kaon propagator Dxiq^) in free space are given, 
respectively, by 

and 

DKiqV) = — Y- 
q\-m\ 

In Eq. (1), o"7rAf^A_ft: should be the isospin-averaged ttA^ KK cross section for an 
off-shell pion, but we replace it by an on-shell one. It can then be obtained from the 
experimental data via 



It has been parameterized by Cugnon et al. and used in Ref . pGf . It can also be calculated 



from theoretical models. We use here the results of Ref. 25 based on a resonance model 



which describes rather well the experimental data. To take into account the off-shell nature 



of the exchanged pion, we introduce as in Ref. p^l ttNAK vertex a pion form factor 

similar to the one used in the nNN vertex, i.e., 

n<i.) - -jT—T^ 

IT J-TT 

with being the cut-off parameter. 

In the kaon exchange contribution, ctkn^kn is the isospin-averaged kaon-nucleon scat- 
tering cross section, and can be determined from experiments by 

I, . 

CkN^KN — -[(^K+p-^K+p + (^K^p^ROp + (^KOp^K+n)^ 

where we have made use of the following relations based on the isospin symmetry: 

O'KOn^ROn = (^K+p^K+p, CFK+n->K+n = CTKOp^ROp, and ax+n^KOp = aK^p^K+n- TheSe CrOSS 



sections are taken from the phase shift analysis of Martin . The off-shell effect is again 
included by the following kaon form factor at both vertices, 

^ yiK) - -j-2 — 
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with Kk being the cut-off parameter. 

For the reaction pp pAK^, where experimental data exist, the cross section is similar 
to Eq. (1) except that the first term on the right hand side is reduced by a factor of 3 due 
to the isospin and axN^KN in the second term replaced by ax+p^K+p- If only the pion 
exchange is included, the isospin-averaged cross section is three times the pp pAK^ cross 
section as in Refs. [i,^]- Since the isospin factor for the kaon exchange is one, the isopsin- 



averaged cross section including both pion and kaon exchanges is thus less than three times 
the pp pAK'^ cross section. 

Next, let us consider the reaction where a kaon is produced in association with a S 
hyperon. Since the S hyperon has three charge states, there are more final states than in 
the case for the A hyperon. The isospin-averaged cross section is given by 

<yNN^NT,K{\fs) = 

Smf^ r{'^^)max 2i f^'^"'^^ J 2fnNN T^ii' 2\ 2|n / 2m2- / \ 

dw.w^Kj dq^-^F (qMDMl a.N^Mw.) 

^ \ / min \^7v / TV 

3m^ms , 2 , /■^^^^+ , 2 /Itve ^.4/ 2 a 2 



i 22 l( ."^^^ dwKwj^kK f J" dq'j^^^F\qj^)q]^\DK{qK)\'^(TKN->KN{wK)- (6) 
Ztt p s J[ujK>^i„ Hik)- rrij^ 

The notations are similar to those in Eq. (1). The mass of the S hyperon and the KNT, 
coupling constant are denoted by ms and fxNT., respectively. The isospin-averaged cross 
section for the 7iN — > T,K reaction is denoted by (Jt^n^j^k and is calculated from 

1 



For the cross sections on the right hand side, we use again the results of Ref. ||2^ based on 
the resonance model which has been shown to fit the experimental data satisfactorily. 

For pp n'E'^K'^, pp pTPK^, and pp pTi^K^, there are experimental data 
available, and their cross sections can be obtained from Eq. (6) by replacing 3(T7r7v^s_ft: in 
the first term by 2a^+p^s+/!'+, ^■kOp-^t,ok+, and o"7r0p^s+_ft'0, respectively, and 3aKN^KN in 
the second term by 2(j^()p^;^+„, (rx+p^K+p, and laROp^K^^p^ respectively. 
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B. iVA ^ NYK 



As mentioned in the Introduction, deltas play an important role in subthreshold kaon 
production from heavy-ion collisions. Unfortunately, there are no experimental data for kaon 
production cross sections from the nucleon-delta and delta-delta interactions. A theoretical 
model is therefore needed to determine these cross sections. In this and next subsections 
we extend the pion and kaon exchange model to calculate kaon production cross sections in 
these reactions. 

For the A^A NYK reaction, the relevant diagrams for pion and kaon exchanges are 
shown in Fig. 2. In each case, there are two different diagrams; one with a kaon produced 
from the meson-nucleon interaction and the other from the meson-delta interaction. For 
A^A NAK reaction, diagram (d) does not exist as there is no A'AA coupling due to 
isospin conservation. The isospin-averaged cross section for this reaction is 

<yNA~*NAK{Vs) = 

t^^,^'^ dw.wiK /tf ^^^%^'(^^)^(^^)l^.(^:)r-.iv^A^(^.) 

m]^m\ ri'^ii)^^^ r(iK)+ 2 Ikna 774/ 2 \ 2 1 n / 2 \|2- / \ /^7^ 

+ .22 ( \ dwKW^kK dq^ — —F [q^)qK\DK[qK)\ crKA^KN^WK). (7) 

The notations are again similar to those in Eq. (1). In the first term, corresponding to Fig. 
2(a), k'^ is calculated from Eq. (2) with rriN replaced by the mass of the delta resonance, m^. 
o'ttA^ak is the isospin-averaged cross section for the vrA —>■ AK reaction which is determined 
from 



1 

— ( 

3 



O'ttA-^AK — -<^n-A++^AK+, 



with <JTr-A++-^AK+ taken from Ref. p5 



In the second term, corresponding to Fig. 2(b), firNA is the vrA^A coupling constant, and 
{qiy± are given by 

{qly_ = m% + ml - 2{ml + p^f/\m], + p'^f'^ + 2pp\ 



{ql)'^ = ml + ml- 2{m\ + v^f'''{rn\ + p'2)i/2 _ 2pp'. (8) 
A(g^) is from the ttA^A vertex and is given by 

In the third term, corresponding to Fig. 2(c), k'j^ is calculated from Eq. (2) with m^r 
and itIt, replaced by ttia and uik, respectively. (Tka^kn is the isospin-averaged cross section 
for the KA KN reaction. Since there is neither experimental nor theoretical information 
on this cross section, we assume that it is the same as the isospin-averaged cross section 
(Jkn^kn for the KN KN reaction. 

Since a delta can decay into a physical pion and a nucleon, the exchanged pion in Fig. 
2(b) can be on shell, which leads to a singularity in the pion propagator. However, a pion in 
nuclear matter acquires a (complex) self-energy due to the strong pion-nucleon interaction. 
The finite imaginary part of the pion self-energy then makes the contribution from the on- 
shell pion finite. Since we are interested in kaon production in nuclear medium, we replace 
the free space pion propagator by an in-medium one, 

Diql) = V-Ff ' (9) 

q^- -U 

with n being the pion self-energy to be specified later. 

The AA^ —>■ NYK reaction includes contributions from both on-shell and off-shell pions. 
In transport models, the on-shell pion contribution is usually treated as a two step process, 
i.e., a delta decaying into a physical pion and a nucleon, and the subsequent production of a 



kaon from the pion-nucleon interaction [^. In this case, we should include only the off-shell 
pion contribution in AA^ NYK. Since it is not clear how one can properly separate the 
on-shell and the off-shell contribution in AA^ —>■ NYK, we shall use in the transport model 
the total AA^ —>■ NYK cross section and neglect kaon production from explicit pion-nucleon 
interactions. A similar strategy has been adopted in Ref. |^ where the eta {rj) production 
cross section from the nucleon-delta interaction has been considered. 
For A^A — s> NTjK, the isospin-averaged cross section is given by 
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<^na^ni:k{Vs) = 

3m% fM^a. , 2,/ f^'''^+ , 2fn 



6mNm^ fK^K)^^ , 2 /™p4/ 2x 2 in ^2x12- / n 

+ A 2 2 ( \ dwRW^k^ dqj^ —F [q^)qK\DK[qK)\ (^ka^kn{,Wk) 

^^^^^^^ /tt''^^x%i^^(g|)^(gi)l^/c(gi)r-i.^^K^(«^i.). (10) 
^TT p s J[wk)^^„ J{ik)'- rriR 

Similar notations as those in Eq. (11) are used. A{qj^), obtained from the KAH vertex, is 

given by 

Mil) = ,0 ^ 3 llK - i^A + mj,f][{ml - m| + q],f - Anilql]. 

The isospin-averaged cross section o't./^^y^k is determined from 

1,4 5 , 



The cross sections on the right hand side are again taken from Ref. p5[. For A^A NT^K, 
a kaon can also be produced from the the KN vertex (Fig. 2(d)). This contribution is given 
by the last term in Eq. (10), where fx as is the KAT, coupling constant, and (gl-)'-]- are 
calculated from Eq. (5) with uin- and tjia replaced by mA and ms, respectively. 

C. AA ^ NYK 

The pion and kaon exchange diagrams for this process are shown in Fig. 3. In the case a 
kaon is produced in association with a A hyperon, the kaon exchange is not allowed because 
of isospin conservation. The isospin-averaged cross section is given by 

^ /'"f ""^ dw.wlK dqlf-i^F\ql)A{ql)\DMl)?-^^^^^ (11) 

On the other hand, both pion and kaon exchanges contribute to the AA — » NJ^K 
reaction. The isospin-averaged cross section is 
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•J [wn ]__ ,__ J{qi) rn^ 



r' -^ dwKwj^k'^ dql^-^F\ql)A{ql)\DK{ql)?^^^^ (12) 

-/[wk] -lipiA lilt,- 



D. The pion self-energy in nuclear matter 

As mentioned earlier, to calculate the pion exchange diagrams for kaon production in 
nucleon-delta and delta-delta interactions requires the inclusion of a (complex) pion self- 
energy. Following Refs. [^ , ^^ -|35| , we calculate the self-energy of a pion with energy u and 



momentum k in nuclear matter by taking into account the delta-hole polarization in the 
random-phase approximation, i.e., 

n(.,k) = ^!^^^^, (13) 

where g' ~ 0.6 is the Migdal parameter due to the short-range repulsion. The pion suscep- 
tibility X is given by 



X(c.,k) ^ ^(^)^^^exp(-2kV5>, (14) 
where p is the nuclear matter density and h ^ Tm-Tr is the width of the form factor [BH 



Including the delta width Fa, oja is approximately given by 

k^ i 

uja ~ 7; h rriA - rriN - -Fa- 

2mA 2 

The delta width in nuclear matter can be determined by extending the delta-hole model to 
include effects of the in- medium pion dispersion relation, i.e., 

TAiu:) =-2/(g,(^)'exp(-2k'V&2) 
■ Im 



j^D^{u,k')j^ + g'^^\, (15) 
where D.^{uj, k) is the in-medium pion propagator given by Eq. (9). 
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The delta decay width and the pion self-energy are then calculated self-consistently. We 
show in Fig. 4 the real (left panel) and imaginary (right panel) parts of the pion self-energy 
as functions of its momentum for a number of energies. The nuclear matter density is taken 
to be 2po, with po = 0.17 fm~^. We see that with increasing pion energy, the real part 
decreases, while the imaginary part first increases and then decreases. 

The pion self-energy is also affected by temperature. For heavy-ion collisions at 1 
GeV/nucleon, the temperature reached in the collisions is below 100 MeV, and accord- 
ing to Ref. its effect is not appreciable. Thus, we shall not consider the temperature 
effect in present study. 

III. KAON PRODUCTION IN BARYON-BARYON INTERACTIONS: RESULTS 

AND DISCUSSIONS 

Our model for kaon production in baryon-baryon interactions involves five coupling con- 
stants and two cut-off parameters. We use the following values for these parameters. 



fnNN 



0.08, 



fnNA 

An 



0.37, 



An 



0.97, 



An 



0.07, 



An 



0.23, 



= 1.2 GeV, A 



K 



0.9 GeV. 



The nNN coupling constant and the pion cut-off mass are typical values used in the Bonn 
model for the nucleon-nucleon interaction |3^. The TriVA coupling constant are taken from 
Ref. [|38| and is consistent with that determined from the A decay width. The KNA and 
KNU coupling constants are the same as those used in Ref. and are close to the 

values used in the Nijmegen as well as the Bonn-Jiirich models for the hyperon- 
nucleon interaction. The KAT, coupling constant is taken from Ref. 0]. Finally, the kaon 



cut-off mass is similar to that used in Ref. p^. The two cut-off parameters have been 
slightly adjusted in order to obtain good agreements with available experimental data on 
kaon production. 
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The comparison of our model prediction with the experimental data for pp pAK~^ is 
shown in Fig. 5, where the solid curve gives the result including both the pion and kaon 
exchange contributions, while the dash-dotted curve gives only the pion-exchange contribu- 
tion. We find as in Ref. |^ that once the off-shell effect at the ttNAK vertex is included, 
the pion exchange alone underestimates appreciably the experimental data. Thus, the kaon 
exchange mechanism is essential for a correct account of the data. As already pointed out 
in Refs. ||T3| , p!6[| , the fact that the pion exchange alone also fits the experimental data in Ref. 



26| is attributed to the neglect of off-shell effects at the ttNAK vertex. 

In Fig. 5, we also show the predictions from the linear parameterization of Ref. 
(dashed curve) and the quartic parameterization of Ref. (dotted curve), which are given, 
respectively, by 

r)A Pmax , 

^linear A'""; 

rriK 

( Pmax \ 4 



with 

P 



CTquartic = 800 ( ) fib, 

TTlK 



1 1/2 

(s - (mN + mA + rnKY){s - {rriN + rriA - itikY) . (16) 



2v^ 

At low energies (below the lowest experimental point of ^/s ~ 2.7 GeV) our prediction is close 
to the quartic parameterization, and is considerably smaller than the linear parameterization. 
Therefore, the use of the linear parameterization for subthreshold kaon production in heavy- 
ion collisions probably overestimates the contribution from the nucleon-nucleon interaction. 
On the other hand, the quartic parameterization greatly overestimates the cross section for 
^/s>2.7 GeV. 

Figs. 6-8 show comparisons of our results with the experimental data from reactions in 
which the kaon is produced in association with a S hyperon. In this case, the pion exchange 
alone can explain the data as the kaon exchange contribution is very small due to a small 
KNTj coupling. This is consistent with the observations in Ref. |^ . 

In transport models, isospin-averaged cross sections are usually used. We denote, by 
<7bb~*nyk, the isospin-averaged kaon production cross section in baryon-baryon interactions 
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including both the NAK and NT^K final states, i.e., (Jbb^nyk = o'bb^nak + o^bb^njik- 
In Fig. 9, we compare aNM->NYK based on the present model (Eqs. (1) and (6)) with the 
parameterization of Ref. (dashed curve) 

a^j, = 72 ^ + 72 ^ (17) 

•rriK rriK 

where p'^^x calculated from Eq. (16) with rri/^ replaced by m^. 

The theoretical results for aNN^NYK are shown for nuclear matter at densities p = 0, po 
and 2po- In general, our results are smaller than the linear parameterization of Ref. P|, 
especially near the threshold. The sudden change of the kaon cross section at ^/s ~ 2.63GeV 
in the linear parameterization is due to the onset of the A^A^ NUK channel. With 
increasing density, the kaon production cross section is seen to increase. This is consistent 
with the findings of Ref. |2^, and is mainly due to the softening of the pion dispersion 



relation, which tends to enhance the cross section. 

Our results for the kaon production cross section in the A^A interaction is shown in 
Fig. 10. The result from the scaling ansatz of Ref. ^ (i.e., three-fourth of Eq. (17)) 
is also shown in the figure by the dashed curve. Our results near the threshold are close 
to the scaling ansatz. The density dependence of o-jq/^^^^YK is, however, different from 
that of aNN^NYK- While the latter increases with density because of the softening of the 
pion dispersion relation, the former decreases with density. This is largely due to the strong 
suppression of the singular on-shell pion contribution at high densities. A similar observation 



has been found in Ref. where the eta production cross section in the nucleon-delta 
interaction is seen to also decrease with increasing density. 

Finally, our results for the kaon production cross section in the delta-delta interaction 
are shown in Fig. 11. Near the threshold they are smaller than the scaling ansatz of Ref. 
(i.e., half of Eq. (17)). Also, the kaon production cross section in the delta-delta interaction 
is considerably smaller than that in the nucleon-delta interaction as a result of the smaller 
vr A YK cross section than the vrA^ YK cross section . Again, we see that the kaon 
production cross section decreases with increasing density. 
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IV. APPLICATIONS TO AU+AU COLLISIONS AT 1 GEV/NUCLEON 



In this section, we shall use the theoretical kaon production cross sections from baryon- 
baryon interactions in a transport model to study kaon production in Au+Au collisions at 



1 GeV/nucleon. This reaction has been previously studied by many groups [p^ -[T7| based 
on various transport models and using mainly the linear parameterization of Ref. for the 
elementary kaon production cross sections in baryon-baryon interactions. 

Ideally, we would like to carry out the calculation in the relativistic transport model as in 
our previous studies |T^|T^. In this way, we can treat properly the medium modification of 
hadron properties |jl3|. Since the kaon production cross sections from baryon-baryon inter- 
actions determined in the present work have been obtained without these medium effects, we 
shall use instead the non-relativistic transport model in which hadron masses are taken to 
be free ones. The difference between the relativistic and the non-relativistic transport model 



for kaon production is, however, insignificant, as already discussed in Ref. [|T7|- The reason 
is simple: in the reaction BB — > NYK, both scalar and vector potentials in the intial and 
final states cancel in leading order. The kaon production threshold in the medium is thus 
not affected by medium effects, which are neglected in the non-relativistic transport model. 
The situation is certainly different for BB — > NNKK and BB NNpp where a net scalar 
potential is present in the final state. Including the change of the production threshold in 
the medium, it has been shown that the relativistic transport model leads to almost two 



order-of-magnitude enhancement of the antiproton yield |^2| and a factor of 4 enhancement 
of the antikaon yield |^3| as compared to those from the non-relativistic transport model. 
In the future, we plan to carry out a more complete study of subthreshold kaon production 
in heavy ion collisions based on the relativistic transport model that includes consistently 
medium effects on the elementary kaon production cross sections and the collision dynamics. 

With a soft Skyrme equation of state, corresponding to a compressibility K = 200 
MeV at normal nuclear matter density, we show in Fig. 12 the energy distribution of 
baryon-baryon interactions that lead to kaon production in central Au+Au collisions at 1 
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GeV/nucleon. There are about 112 baryon-baryon collisions with available energies above 
the respective kaon production thresholds. The average energy in these collisions is about 
2.64 GeV. Since the threshold energy for the BB AYK reaction with a delta in the 
final state is about 2.84 GeV, our neglect of its contribution to the kaon yield is therefore 
justified. We would like to point out that an average energy of 2.64 GeV corresponds to a 
kaon maximum momentum of 0.275 GeV, which is very close to the average kaon maximum 
momentum of 0.272 GeV obtained in Ref. [jl^] based on the relativistic transport model 



and a soft (relativistic) equation of state. This is an indication of the similarity between the 
relativistic and non-relativistic transport model descriptions of kaon production in heavy-ion 
collisions. 

The kaon production probabilities from different reactions are shown in Fig. 13. The 
results based on the theoretical elementary cross sections are shown in the left panel. The 
contribution from the nucleon-delta interaction is most important and accounts for about 
65%. As in previous studies, deltas play the most important role as reactions involving delta 
resonances produce about 90% of the kaons. 

We have also carried out a VUU calculation using the Randrup-Ko parameterization for 
kaon production in baryon-baryon interactions. The results are shown in the middle panel 
of Fig. 13. We see that the kaon yield in this case is only about 20% less than that obtained 
with the theoretical cross sections calculated in the present paper. 

To see quantitatively the difference (or similarity) between the relativistic and non- 
relativistic descriptions of kaon production in heavy-ion collisions, we reproduce the Fig. 
8 of Ref. |]13[ in the right panel of Fig. 12. The relativistic transport model is seen to 



lead to a reduction of about 30% in the total kaon yield. This reduction is mainly due to 
the momentum-dependent nuclear mean-field potential included in the relativistic transport 
model. This momentum- dependent potential gives rise to an additional repulsion besides 
that from the compressional pressure, and a smaller maximum density is thus reached in the 
relativistic transport model. This reduction in kaon yield is much smaller than that found 
in Ref. [§] based on the non-relativistic QMD model with a momentum-dependent Skyrme- 
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type potential, where the reduction factor is about 3. The reason for this difference between 
these two calculations has been discussed in Ref. |T^. We believe that the treatment of the 
momentum-dependent mean-field potential in Ref. |^ is incomplete as it has not taken into 
account the difference in the intial and final potential energies in the BB —>■ NYK reaction. 
Since baryons in the initial state have larger momenta, the intial potential energy in the 
above reaction is larger than that of the final state and can thus be used for kaon production. 
If this potential difference is properly included, as in the relativistic transport model, the 
difference between the kaon yields obtained in QMD calculations with and without the 
momentum-dependent potential should be small. 

Finally, we show in Fig. 14 the kaon momentum spectrum obtained in this study. The 
theoretical result is seen to overestimate the experimental data |^ by about 30%. As men- 
tioned earlier, the relativistic transport model would lead to a reduction of a similar mag- 
nitude in the kaon yield as compared to the non-relativistic one. We thus expect that a 
consistent calculation based on the relativistic transport model and using the theoretical el- 
ementary cross sections which properly include the medium effects, will give a better account 
of the experimental data. 



V. SUMMARY 

In this work, we have extended our previous studies of kaon production in Au+Au 
collisions at 1 GeV/nucleon by using in the transport model the elementary kaon production 
cross sections calculated from a theoretical model. The main purpose of this study is to 
examine our previous conclusions concerning the nuclear equation of state at high densities 
| T7| and the kaon scalar potential [|13| in dense matter. A particular concern is the possible 



underestimate of the kaon production cross section from the nucleon-delta interaction in 
the scaling ansatz of Randrup and Ko. For this purpose, we have constructed a one-pion 
plus one-kaon exchange model for kaon production in baryon-baryon interactions. The 
parameters of the model are fitted to available experimental data in the proton-proton 
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interaction. This model is then extended to the nucleon-delta and delta-delta interactions. 
We have found that near the kaon production threshold, which are relevant for subthreshold 
kaon production, the cross sections in the Randrup-Ko parameterization and scaling ansatz 
that have been used in previous studies by us [0,113] s-nd by other groups [0,|1^,0 are 



overestimated for the nucleon-nucleon interaction but underestimated for the nucleon-delta 
interaction. 

We have then used these theoretical cross sections to calculate in a non-relativistic trans- 
port model the kaon yield in Au+Au collisions at 1 GeV/nucleon. The results are found 
in reasonable agreement with experimental data. We expect that a consistent relativistic 
transport model calculation will lead to an even better agreement with the data. The con- 
clusion of this extensive study is thus clear: we need both an attractive scalar potential 
and a soft nuclear equation of state to account for the experimental data from the Kaos 
collaboration at the GSI. 

We are grateful to J. Aichelin, G. E. Brown, B. A. Li, and K. Tsushima for helpful 
discussions. This work was supported in part by NSF Grant No. PHY-9212209 and the 
Welch Foundation Grant No. A- 11 10. 
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Figure Captions 

Fig. 1: One-pion and one-kaon exchange diagrams for tlie A^A^ — > NYK reaction. 
Fig. 2: Same as Fig. 1, for the A^A — NYK reaction. 
Fig. 3: Same as Fig. 1, for the A A NYK reaction. 

Fig. 4: Real (left panel) and imaginary (right panel) parts of the pion self-energy in nuclear 
matter with density p = 2po- 

Fig. 5: Comparisons of model predictions with experimental data for the kaon production 
cross section in pp pAK^. The dashed-dotted curve gives the one-pion exchange 
contribution, while the solid curve includes both the one-pion and one-kaon exchange 
contributions. The dashed curve is based on the linear parameterization of Randrup 
and Ko [Q, while the dotted curve is based on the quartic parameterization of Zwer- 
mann B . The experimental data are from Ref . . 



Fig. 6: Same as Fig. 4, for pp pTPK'^ . 
Fig. 7: Same as Fig. 4, for pp pTj^K^. 
Fig. 8: Same as Fig. 4, for pp nJ^^K^. 

Fig. 9: Isospin-averaged cross section aNN-^NYK{= o-^n^nak + crNN^NY:K)- The solid 
curve is our model prediction and the dashed curve is from the parameterization of 
Ref. i. 

Fig. 10: Same as Fig. 9 for aNA^NYK{= (^na^nak + (^na^nek)- 
Fig. 11: Same as Fig. 9 for (Taa^nyk{= (^aa^nak + <yAA~.NSK)- 

Fig. 12: The energy distribution of baryon-baryon collisions that are above the kaon 
production threshold. 
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Fig. 13: Kaon production probabilities from different channels obtained in three calcula- 
tions. 

Fig. 14: Comparisons of kaon momentum spectra obtained in this work using the theo- 
retical kaon production cross sections (solid curve) with the experimental data (open 
squares) from Ref. [Q. 
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